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5t , We study string dynamics in the early universe. Our motivation is the pro- 

posal of Brandenberger and Vafa, that string winding modes may play a key 
role in decompactifying three spatial dimensions. We model the universe as a 
homogeneous but anisotropic 9-torus filled with a gas of excited strings. We 
adopt initial conditions which fix the dilaton and the volume of the torus, 
but otherwise assume all states are equally likely. We study the evolution 
of the system both analytically and numerically to determine the late-time 
behavior. We find that, although dynamical evolution can indeed lead to 
three large spatial dimensions, such an outcome is not statistically favored. 



1 Introduction 



An enduring challenge for string/M-theory is to provide a more complete 
picture of the early universe than has been found using conventional, point- 
particle approaches. To this end, a growing body of research has studied the 
dynamics of strings and branes in a cosmological setting, as opposed to the 
more widely investigated case of a static background. Intriguing but as yet 
incomplete results have been found for higher dimensional cosmologies and 
cosmologies based on braneworlds. These include mechanisms for resolving 
or avoiding cosmological singularities, and for generating subtle modifications 
to the primordial microwave background power spectrum. Further progress 
on these key theoretical and observational issues, however, requires a more 
refined grasp of the dynamical properties of strings and branes when sub- 
ject to extremes of temperature, density, and curvature. The current paper 
provides a modest step in this direction. 

The formalism we develop can, in principle, be applied to a wide range 
of string cosmology questions. But following our earlier works PP our 
immediate goal is to find a dynamical mechanism within string/M-theory 
that generically gives rise to a universe with precisely three large spatial 
dimensions, with all other spatial dimensions unobservably small. Such an 
asymmetric dynamical evolution is perhaps the most basic task of string/M 
cosmology. However, a decade and a half after the first attempt, no satisfac- 
tory picture has yet emerged. 

By way of brief history, in [3] and 0] the authors made use of T-duality 
in a (spatially) toroidal universe to argue that strings wound around nontriv- 
ial cycles impede the growth rate of the spatial dimensions they wrap. The 
fastest expansion will therefore be achieved by dimensions that shed all their 
winding modes through string winding/anti- winding annihilations. Because 
string worldsheets are two-dimensional, pairs of strings will generically inter- 
sect in four or fewer spacetime dimensions, leading 3j and ^ to argue that 
at most three spatial dimensions will shed their windings and subsequently 
grow with time. Various aspects of this proposal have been investigated and 
generalized. In [Sj, a numerical study of a gas of strings in a toroidal uni- 
verse was carried out and the naive dimension counting argument for string 
annihilations used in [3j (2+2 = 3+1) was verified in a static background. 
In jS], the analysis was further extended to simply connected toroidal orb- 
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ifolds, and it was argued that pseudo-winding modes with sufficiently long 
lifetimes could allow the arguments of [3| to apply in phenomenologically 
relevant backgrounds. In fTI , higher dimensional branes were included in the 
analysis, without invalidating the conclusions of jS] . Other studies of string 
windings and brane gases can be found in |5]-|24j. 

All these works, however, fail to account for the detailed cosmological 
dynamics. In [T] we partially addressed this deficiency by using eleven- 
dimensional supergravity to study cosmological evolution in the presence of a 
brane gas and found encouraging results: for a suitable configuration of brane 
wrappings suggested in JT] and |1| and based on naive dimension counting 
arguments, the dynamics does indeed drive an asymmetric evolution yield- 
ing three large dimensions. In [2^ we went further and studied the coupled 
Einstein-Boltzmann equations for a thermal brane gas and found that despite 
the naive dimension counting arguments, only highly specialized initial condi- 
tions yield the desired brane wrapping configuration. In particular, we found 
that the spatial expansion driven by the brane gas is generically too fast for 
brane interactions to generate the expected anisotropics; instead, the branes 
quickly freeze out. However, this analysis still held out the possibility of a 
loop-hole that would allow one to evade this discouraging conclusion. In the 
string theory corner of M-theory moduli space (the very scenario studied in 
Brandenberger and Vafa's initial paper the growth of spatial dimensions 
appeared to be slower, perhaps allowing sufficient string-string interactions 
to yield the desired asymmetric winding configuration and hence asymmet- 
ric expansion. The main purpose of this paper is to study this possibility in 
detail. 

In section 2 we set up the basic framework of dilaton gravity, and in sec- 
tion 3 we discuss the equilibrium thermodynamics of a string gas. In section 4 
we introduce the Boltzmann equations which govern string annihilation and 
give a preliminary discussion of the phenomenon of freeze-out. In section 5 
we describe our method for sampling from the possible initial states of the 
universe, and discuss holographic bounds on the space of initial conditions. 
In section 6 we present numerical results, and show that other than for a 
narrow range of initial conditions we get "all" or "nothing" evolution: either 
there are too few strings to keep any spatial dimensions small, or string inter- 
actions freeze out so quickly that a large number of wrapped strings survive 
and prevent any dimensions from growing large. 
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2 Dilaton gravity 



We start with type II string theory compactified on a 9-torus, with metric 

9 

ds^ = -de + a'Y, e^^^^'^de^ < < 27r . (1) 

i=l 

The action for dilaton gravity is 

5 = -V / d^'x^e-'^ [R + mr + ■■■) (2) 

where (Polchinski |25|, 13.3.24) Kj^ = |(27r)^(a')^ From now on we set 
a' = 1. Following Tseytlin & Vafa |5, we define the shifted dilaton 

ip = 2<P-J2>^^ (3) 

i 

SO that the action reads 

S={2nr [dte~-{j2>^"-^")- (4) 

When one couples dilaton gravity to a matter system the time-time compo- 
nent of the Einstein equations yields the Hamiltonian constraint (or Friedman 
equation) 

{2nfe--{^^-J2)^)=E (5) 

i 

where E is the total matter energy. This constraint implies that 0^ never 
vanishes; we choose the direction of time so that < 0. The dilaton equation 
of motion is 

i 

The scale factors obey 

A, - = ^e^Pi (7) 

where the "total pressure" Pj = — 1£- is obtained by varying the matter free 
energy with respect to Aj. Pi is equal to the ordinary pressure in the i*^ 
direction times the spatial volume. 
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Note that T-duality in the ith direction takes the simple form 



Aj — s> — Aj, (f invariant . 

This leaves the dynamics unchanged, provided that E is invariant and Pj 
changes sign. 

To get oriented we first study the vacuum equations, with all pressures 
set to zero. The equations of motion reduce to 

i 

Xi-ipX, = 0. (9) 

Besides the trivial solutions in which the dilaton and radii are constant, a 
Kasner-like branch of solutions can be obtained as follows. If the pressures 
vanish the energy E is conserved. Then (jH} can be reduced to an equation 
just for if, with general solution 



(fit) = log 



16tt^/E 
t{t + C) 



(10) 



(we have suppressed one constant of integration corresponding to an arbitrary 
shift in t). One can then integrate the Aj equations of motion to find 

Ut) = A + Bdogj^. (11) 

The constants of integration Ai are arbitrary, while in order to satisfy the 
Hamiltonian constraint Bi and C must satisfy (1 — Bf) = 0. Thus 
either C = and the radii are static, or V. = 1 and the radii are time 
dependent. In both cases, the dilaton rolls monotonically towards weak cou- 
pling. 

We now turn to the the late-time asymptotic behavior of solutions to the 
dilaton-gravity equations.^ First suppose the pressure is negligible. Pi ~ 0, 
as is the case for a universe in equilibrium with all radii sufficiently close to 
the self-dual radius. At late times the universe will approach the Kasner-like 
solution (fTUj). (fTTj). with the asymptotic behavior 

const. X . 

e'P ~ — — e^' ~ const. (12) 



similar analysis was performed for M-theory in 
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Thus if Pi ^ the dilaton rolls monotonically while the radii approach 
constants in string frame. ^ 



Now suppose that at late times we have m unwrapped dimensions ■ ■ ■ 
and 9 — m wrapped dimensions x™""^^ ■ ■ - x^.^ Without loss of generality we 
can go to a T-dual frame where the unwrapped dimensions are all larger 
than string scale. In this frame we expect that at late times the universe 
will be dominated by a radiation gas in the unwrapped dimensions; the pres- 
sures should vanish in the wrapped dimensions due to a cancellation between 
winding and KK modes. That is, we expect 



An ansatz which captures the late-time behavior is 
1 .x.. it' ^ ^ 



const. t = m + 1, . . . ,9 

Plugging this ansatz into the equations of motion fixes 

2m 2 

« = ^ P = ^- (15) 

m + 1 m + 1 

The dilaton rolls monotonically to weak coupling, while the unwrapped di- 
mensions grow with time and the wrapped dimensions have fixed sizes. Thus 
if the string winding dynamics in the early universe favors m = 3, as sug- 
gested by the dimension counting argument reviewed in the introduction, one 
could naturally explain why three spatial dimensions become large. 



3 Equilibrium thermodynamics 

In a coupled matter /gravity system the matter energy is determined by the 
Hamiltonian constraint (0). For stringy matter Tseytlin & Vafa long ago 
presented a simple picture of the corresponding thermodynamics which is 
suitable for our purposes. There are two possible phases. 

^This was also shown in section 5 of PP. To relate the two solutions note that PP worked 
in terms of M-theory time Im, related to the string- frame time used here by ts ~ ^m^- 

■^A dimension is called unwrapped if > and the number of winding strings 
vanishes, or if Aj < and the number of momentum modes vanishes. 
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3.1 Hagedorn phase 



Strings have a limiting Hagedorn temperature [22] • For weakly-coupled type 
II strings the limiting temperature is Th = Near this temperature the 
canonical ensemble fails and one must use the microcanonical ensemble |2Zj • 
In the Hagedorn phase the universe contains a dense gas of winding and KK 
modes. To a good approximation the free energy F = E — ThS vanishes, so 
the microcanonical entropy is given by S{E) = E/Th- The total pressure 
also vanishes, Pj = — 1£- = 0. 

The thermodynamics of strings in the Hagedorn phase has been studied 
by Deo, Jain and Tan |2H1 121] • They employ the microcanonical ensemble, 
with a fixed energy E and a fixed net winding charge in the universe; in a 
compact space the latter vanishes. They show that the average number of 
type II strings present with winding charge vector w and energy e is given 
by 

D{e, w, E) = -u{e, ^)rf/2e-«(^.^)w^^-'w/4 ^^g^ 



where 



e 



VdetA 

u{e,E) 



e{E - e) 
A - ^ 5 

Here Ri = e^^ . As a consistency check, note that the total amount of energy 
in strings indeed adds up to E: 



de j d'^w eD{e,w,E) = E. 



We will ignore diagonally-wound topologies (where a string is simultane- 
ously wound on several dimensions) and assume that we have 9 unidimen- 
sional string gases. That is, we set o? = 1 and assume that the total energy 
available to each dimension is 1/9 of the total energy in the universe. Thus 
the distribution for a single winding charge Wi is given by 

D{e,Wi,E) = ^^v^M(e,E/9)exp [-u(e, E/9)u7- vr^P. ] . (17) 
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The (thermally averaged) total number of positive windings Wi is then given 
by 

{W,) = de dw, WiD{e,w,,E) = ^ . (18) 
Jo Jo Vl^-nRi 

Note that in (fTHj) we are only counting positive winding in some direction; 
the net winding is zero. The total length of string present is proportional to 
the energy while the physical distribution of string on the torus amounts 
to a random walk. Thus the average positive winding is simply the average 
distance from the origin of a random walk. This goes as the square root of 
the length, or equivalently the square root of the energy. 

We have computed the average number of winding modes, but a similar 
result holds for the average amount of positive KK momentum present, just 
by replacing Ri — > XjRi. 



3.2 Radiation phase 

Below the Hagedorn temperature the string oscillators make a negligible 
contribution to the partition function, so we can focus on single-string states 
which are labeled by an integer-valued momentum vector rii and an integer- 
valued winding vector Wi. In the absence of a i?-field the string energy levels 
are 




The corresponding free energy for a gas of strings is 

pF = 128 logtanh(/5e(n,w)/2) 

nw=0 

where we have taken into account that for the type II string we have 128 
bosonic and 128 fermionic species of excitations. The condition n ■ w = 
enforces level matching. For dimensions which are large compared to a 
thermal wavelength we can approximate momentum sums by integrals and 
neglect winding; likewise for dimensions which are small compared to an 
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inverse thermal wavelength we can approximate winding sums by integrals 
and neglect momentum. Any remaining intermediate-sized dimensions are 
frozen, with no excitations. Thus we have an approximate expression for the 
free energy 

/3F ^ 128 n 27ri?. H ^ / 7^ logtanh(/3|p|/2) (20) 

large small 

where d is the total number of unfrozen dimensions. At this point it is 
convenient to order 

|Ai| > IA2I > ■ ■ ■ > IA9I (21) 
and to define the T-duality invariant 'volume' of the d unfrozen dimensions 

d 

Vd = \{27ie\''^\ . (22) 
1=1 

Then (j20|l can be identified with the free energy of a massless ideal gas in d 
spatial dimensions in a box of volume Vd- To write an equation of state we 
use the fact that in d spatial dimensions an ideal gas has an energy density 
p = CdT'^'^^, where (for 128 bosonic and 128 fermionic degrees of freedom) 

_2dK{d+l)_ , 
The energy, entropy and total pressure of the gas are given by 



E = CdVdT''+' 
+ 1 

IT' 

OF _ ( sign{Xi)E/d i = l,...,d 



d 



dXi 1 i = d+l,...,9 

However we still need to determine d. If the energy is very small then all 
dimensions are frozen. As the energy increases Ai,A2,... will successively 
unfreeze. This means that the temperature in the radiation phase is given 
by 

Trad = min — - . (24) 

\CkVkJ 
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The value of k which minimizes the right hand side is equal to the number 
of unfrozen dimensions. Trad calculated in this way could be larger than 
the Hagedorn temperature; this signals that the system is actually in the 
Hagedorn phase. That is the true temperature of the system is min(Trad, Th). 

In the radiation phase we can compute the amount of positive KK mo- 
mentum present in equilibrium by using the fact that for a one- dimensional 
massless gas {Ni) is related to the pressure by 

(Ni) = ^PiRi. (25) 

This estimate makes sense for Ri ^ V^, in which case we also have (Wi) — 0. 
If on the other hand Ri <C a/o' we just use the T-dual formulas 

(W,) = -\P^/R^ (26) 

m = 0. 



4 Winding and KK Annihilations 
4.1 Boltzmann equations 

If the universe was in equilibrium wc could just insert the results of the previ- 
ous section into the dilaton-gravity equations of motion. But the Brandenberger- 
Vafa scenario is driven by departures from thermal equilibrium. A crude way 
to keep track of these departures is to let Wi be the amount of positive wind- 
ing charge around dimension i. Likewise let Ni be the amount of positive 
Kaluza-Klein momentum in direction i. Of course since the space is com- 
pact there must also be units of anti-winding charge and iVj units of 
anti-KK- momentum. 

Let us, for the time being, assume that each unit of charge is carried by 
a single string: that is, there are strings each wound once with positive 
orientation, similarly for the KK modes. We will also, for the time being, 
assume that the strings have no oscillator excitations. Then the annihilation 
of these momentum and winding modes is governed by Boltzmann equations. 
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similar to the equations that govern the evolution of M2-brane winding j2] : 



m ^ _(M)e^-2A.(iv2_(iv,)2) (27) 
at zn 

^ = _(„,•)=) 
at Zn 

The cross-section for bosonic wound strings was calculated by Polchinski jHO] , 
who found f{v) = 2/(1— f^) for two anti-parallel strings moving with velocity 
V. Subsequent studies have evaluated this quantity for F- and D-strings [HT] . 
We will set f{v) ~ 2, appropriate for a gas of slowly- moving strings. 

Let us make a few comments on the structure of these Boltzmann equa- 
tions. First, note that they are invariant under T-duality. Second, note that 
they respect the dimension-counting arguments of Brandenberger & Vafa 
An implicit factor of the inverse volume of the universe is present in the 
definition of e*^ But due to the factor e'^^^ upstairs in the equation for 
dWi/dt, strings wrapped on a large 3-torus will still be able to annihilate 
effectively, just like particles moving in one large spatial dimension. Finally, 
we should contrast our Boltzmann equations with the results presented in P , 
which were appropriate for strings with a dilaton-independent cross-section 
such as cosmic strings. 

Continuing our study of the Boltzmann equations, we now include string 
oscillator excitations but still restrict attention to unit winding and mo- 
mentum charges. With oscillators excited a more accurate cross-section is 
obtained by replacing^ 

exp(— Aj) e,, 

in the Boltzmann equation for dNi/dt, where Cj is the average energy of a 
string with a unit of KK momentum. Likewise in the dWi/dt equation we 
should replace 

exp(-FAi) Si 

where Si is the average energy of a unit winding string. This is supported 
by the results of Lizzi and Senda [221 , who redo Polchinski's calculation for 
two highly excited strings - which have many oscillator excitations but no 
winding - and show that the interaction rate goes like the product of the 

^We thank R. Myers for bringing this issue to our attention. 
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energies of the two strings. This modifies the Boltzmann equations to read 

at 71 

^ = --e-{6.nW^-{Wn (28) 
at TT 

We may estimate the typical energy per momentum or winding mode from 
the distributions given in the previous section: 

E 

Hagedorn phase : (cj) = tttttt (momentum modes) 

E 

('^») = n/Ti/\ (winding modes) (29) 

Radiation phase : (e^) = (momentum modes) 

{6i) = Ri (winding modes) . 

Here we have assumed that in the Hagedorn phase the energy is equally dis- 
tributed between dimensions. Note that in the Hagedorn phase the average 
energy per mode scales as a/E. 

Finally we consider strings that are multiply wound around each dimen- 
sion, so that the number of positively-wound strings and the winding charge 
Wi are not necessarily the same. We may think of the winding charge as 
made up of Wi open unit strands that are braided together to form closed 
strings. Depending on the braiding, there can be anywhere from 1 to 
closed strings present. Also depending on the braiding, an individual closed 
string can carry anywhere from 1 to Wi units of winding charge. There are 
Wi\ ways of braiding the strands; we assume all braidings are equally likely. 
Then the typical strand is part of a closed string that carries winding charge 
[Wi + l)/2. The cross-section of a string is proportional to its length and 
hence winding charge, thus we expect the typical string-string cross-section 
to be enhanced by a factor {{Wi + l)/2)^. Making this modification to the 
cross-section, and rewriting the Boltzmann equation as an equation for the 
rate of change of the positive winding charge, a net enhancement factor of 
{Wi + l)/2 appears on the right hand side relative to (f2H|) . 

With an analogous modification for multiple-momentum strings the Boltz- 
mann equations read 

dN, {Ni + 1) 



dt 27T 



^{e,Y{Nt-m') (30) 
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dWi 



m + 1) 



dt 



where the average energies per strand are still given by (jSH))- 

One could question our assumption that all braidings are equally likely. 
Although this seems like a reasonable assumption when the universe is small, 
on entropic grounds it could be that as some dimensions grow large singly 
wound strings become favored. Rather than study this issue directly, in our 
numerical work we will investigate the two extreme possibilities: all strings 
singly- wound as in (^Hj). or all braidings equally likely as in (jHUj) . 



4.2 Freeze-out 

In an expanding universe the evolution of a species depends on the species' 
annihilation rate F and the cosmological expansion rate (or Hubble param- 
eter) R/R = A. For dilaton gravity one also needs to take into account the 
rate of change of the dilaton 0. To determine whether a nonzero number of 
strings survive to the asymptotic future, we need to study how these param- 
eters evolve. If the annihilation rate of wound strings decreases too rapidly it 
could undermine the naive dimension counting arguments, which implicitly 
assume that F remains non-zero. 

To illustrate the possibility of freeze-out consider the following simple 
situation. At some initial time set all Aj = 0. Introduce the same number 
of unit winding and unit momentum strings in all directions: Ni = Wi = N. 
Suppose further that no oscillators are excited. Then the pressures Pi all 
vanish, and it is consistent to set the logarithmic scale factors Aj = for all 
time. The remaining equations of motion are very simple. As explained in 
section 2, the shifted dilaton obeys ip = |(^^ with solution 



Here A and to are two constants of integration. As expected the dilaton rolls 
monotonically to weak coupling. The Boltzmann equation (for singly-wound 
strings) reads 



A 
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To get a feel for whether the strings will freeze out it suffices to set (A^) = 0. 
Then the general solution is 

As long as the integral stays finite as t — oo a non-zero fraction of the strings 
will freeze out. In the case at hand the pressure vanishes, which means the 
total energy in matter does not change with time; since the radii are fixed 
the average energy per string (e) also remains constant. Then the integral is 
strongly convergent, and 

1 _ 1 A{e)^ 
N{t) ~ n{t - to) ■ 

Here we have defined 1/A^oo = ^/N{ti) + A{ef/Ti{ti - to). As t ^ oo a 
non-zero fraction of the unit winding and unit momentum strings do freeze 
out, with N{t) Noo. 

This sort of behavior should be fairly generic, even for solutions that do 
not sit precisely at the self-dual radius. The pressure vanishes as long as 
one remains in the Hagedorn phase, giving give rise to a conserved matter 
energy. Moreover, if the radii change slowly with time, the average energy 
per string still remains roughly constant. The dilaton, however, will still 
roll monotonically towards weak coupling, and as long as it does so quickly 
enough for the integral in pip to converge, some strings will freeze out. The 
enhanced cross-section due to multiple winding in (jHUj) does not change this 
outcome. 

This is troubling for the Brandenberger-Vafa scenario, as it shows that 
simple dimension- counting arguments can fail to capture the true dynamics 
of winding strings. Our goal in the remainder of this paper is to undertake 
a detailed numerical investigation of the likelihood of freeze-out. 



5 Initial conditions and Holography 

We would like to choose initial conditions at random, so as to uniformly 
sample the possible states of the early universe. In practice we proceed by 
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fixing the initial value of the shifted dilaton and the initial volume of the 
universe V. All other degrees of freedom will be given random initial values, 
drawn from the probability distribution worked out below. 

Ideally, we would average over all possible values of the "coordinates" Aj, ip 
together with their canonical momenta using the Liouville measure obtained 
from the action The microcanonical volume of phase space for dilaton 
gravity plus matter is 

n J d^Xd^Xdifidifi e-^^^e^ (32) 

where 5* is the matter entropy. In the Hagedorn phase this is given by 

S = E/Th = {2nye-^{^' - ^ X',)/Th . (33) 

i 

Thus at the level of supergravity the initial conditions which maximize the 
entropy are 

(p —>■ —oo weak string coupling and large volume 

^ — oo effective coupling rapidly decreasing (34) 

Aj = constant size of torus 

To set initial conditions we first fix a value of if. In order for effective su- 
pergravity to be valid we must have e*^ -C 1. Note that since we're working 
with effective supergravity, not string theory, only the value of the shifted 
dilaton matters and we don't need to worry about the underlying dilaton (j) 
defined in Q becoming large. Next we fix a value for ip. For supergravity to 
be valid we must have > — 1. Then from (|33|) note that the A are Gaussian 
distributed, with a characteristic spread 

(AA,)2 = T^,(27r)-V. (35) 

For simplicity we take the Aj's to be uniformly distributed about zero, in 
the interval [-v/T^(27r)-^e'^/^ V^(2vr)-^e'^/2]. Note that for < we'll 
have — 1 < Aj < 1. The scale factors Aj do not appear in the entropy, so we 
take them to be uniformly distributed, subject only to the constraint that 
the T-duality invariant 9-volume V defined in (I22|l has the specified value: 
E.|A,,|=log(l^/(2vr)9). 
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Following the M-theory analysis of we can ask about the holographic 
bound ISSl El ESI Ei • This is 



As _ 271 Ae _ 2nnsR% 
^-4G" - i(27r)7(a')4 ^"^^^ 

where fig is the area of a unit and the subscript E reminds us that this 
must be calculated in the Einstein frame. We wish to convert this to string 
frame, with (note that the regular dilaton, not ip, appears below) 

Re = e-^'^Rs (37) 

making 

2fis7r^ 

S < pl^e-- (38) 

on an isotropic torus. Comparing to the initial value we get a bound (again 
assuming we start in the Hagedorn phase) 

S = nV8{2nre-- (^0^ - < (39) 

^V8{2nW < ^e-^ (40) 

2ne^ < (41) 
~ 128^202 

We can interpret this bound 

V ~ (2™^)» < ( -A-l ' = ( —4-) ' (42) 



or equivalently 



128y20V V42Oy20^ 



. 0.16 , , 

^^VTJ^- (43) 
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6 Numerical Analysis 



6.1 Initial Conditions 

Our simulations proceed by generating multiple sets of initial data, solving 
the equations of motion numerically, and looking at the number of wrapped 
dimensions at late times after freeze-out has taken place. In each run we 
fixed the initial values of V and (p. For the most part we started with 
if = —1; this maximizes the entropy while keeping supergravity valid. ^ All 
other initial conditions are allowed to fluctuate randomly. The A, are chosen 
from the flat distribution described in the previous section. The initial A, 
are generated by choosing nine random numbers in [—1, 1] and applying an 
overall scaling so that the initial volume matches the specified value. To 
assign initial values to Ni and Wi we compute the mean values from section 
3 and then add random thermal fluctuations about the mean, of magnitude 

AN, ^ vW AWi ^ ^/{Wii ■ (44) 

Note that we do not impose the holographic bound ()43|) on our initial data; 
we have some comments on this below. 

To evolve the system we use a Runge-Kutta algorithm. At each time step 
we begin by computing the total matter energy from the Hamiltonian con- 
straint (0). The equilibrium thermodynamics discussed in section 3 enables 
us to decide whether the system is in a Hagedorn or radiation phase. Based 
on this we compute the corresponding thermal expectation values (Ni), (Wi). 
We then use the equations of motion to evolve to the next time step. 

The dilaton-Einstein equations of motion are given in (0). To solve 
them we need an expression for the pressures Pi. We set^ 

= 2 {N,e-^^ - W^e^^) . (45) 

^Note that we are considering the case where the dilaton is roUing towards weak cou- 
pling. 

^Here we are relating the actual pressures Pi to the actual values of A^^ and Wi, allowing 
for departures from thermal equilibrium. Thus H45|l should not be confused with (|25|l . (|26|l 
where we used the equilibrium pressures to compute the thermally averaged values of Ni 
and Wi- Of course for a system in equilibrium the expressions are compatible. 
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This simple estimate is valid when no oscillators are excited; thus it should 
be accurate in the radiation phase. In the Hagedorn phase oscillators are 
excited and the pressure receives corrections. However given the equilibrium 
values of Ni and Wi in the Hagedorn phase (fTSj) . (fT^ note that on average 
a cancellation makes the pressure vanish. Thus, although it would be nice 
to have a more precise expression for the pressure, we do not expect any 
refinements to ()45|) to significantly affect our results. 

The Boltzmann equations were discussed in section 4. To allow for the ef- 
fects of multiply- wound strings we ran simulations using two different versions 
of the Boltzmann equations given in (j^H|) and (jHUj) . The first is appropriate 
for strings that only carry one unit of winding or momentum charge, while 
the second is appropriate for strings with multiple winding or momentum 
charges. 

In Figure 1 we show what happens when we vary the initial values of ip 
and V, starting with the initial condition = — 1. We show the average 
number of wrapped dimensions present both in the initial configuration and 
after freeze-out. Clearly a final state with three unwrapped dimensions is 
not dynamically favored. If one begins at reasonably strong coupling then 
very few strings are present in the initial state, while if one begins at weak 
coupling string interactions turn off too rapidly for the required annihilations 
to occur. In either case three large spatial dimensions is not the most likely 
late-time geometry. 

We have explored what happens if the initial value of ip is decreased, since 
the holographic bound (|43|) restricts the allowed values of this quantity. In 
Figure 2 we show the initial and final number of wrapped dimensions starting 
with the initial condition ip = —0.15. The qualitative outcome is the same, 
just shifted to more negative initial values of (p. This is not surprising, given 
the Hamiltonian constraint (0): roughly speaking a change in 0^ can be 
compensated by shifting ip so as to keep the total energy fixed. 

In Figure 3 we show the distribution of initial winding configurations for 
universes that end up with three large dimensions. Although the number of 
wrapped dimensions can either increase or decrease with time, it is unlikely 
that one can begin deep in the Hagedorn phase with nine wrapped dimensions 
and end up with a three dimensional universe. 

In Figure 4 we show how the distribution of final winding configurations 
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depends on the initial value of ^p. Although for a rather narrow range of 
three dimensions is the favored outcome, the distribution of final dimension- 
ality is not very sharply peaked. 

All results presented so far have been based on the multiply-wound cross 
section (|3(jp. We have studied what happens if we evolve the system using the 
singly- wound cross section ()28|) . The change in the results is negligible, much 
less than the widths of the distributions shown in Fig. 4. Thus the qualitative 
outcome is the same, with no dynamical preference for three dimensions. 



7 Conclusions 

Our results indicate that - within the context of our approximations - the ex- 
pansion of the universe has an "all or nothing" character. If initial conditions 
are such that one begins with many wrapped strings, the strings typically 
freeze out and keep all dimensions small. On the other hand if one begins 
with few wrapped strings, the strings typically annihilate and all dimensions 
decompactify. Between these extremes there are initial conditions that lead 
to three large dimensions, but such initial conditions are not generic. Fine- 
tuning the initial conditions to yield three large dimensions is thus possible, 
but runs counter to the goal of the string gas program: finding a mechanism 
in which generic initial data yields three large spatial dimensions. 

The unexpected chink in the Brandenberger-Vafa scenario we have found 
is that due to the rolling dilaton the string annihilation cross section be- 
comes weaker than previously realised. To avoid this impasse we would need 
a mechanism for keeping the string annihilation cross section sufficiently ro- 
bust.'^ There are ways in which this might be accomplished (e.g. strings 
wound on finite fundamental groups ||6j, in confining backgrounds ^Slj, or 
having unusual kinematic configurations |SZ1), but as yet none have been 
studied in adequate detail to determine their viability. Also we should note 
that, even if one manages to stabilize the cross- sect ion, one would still have 
to face the issue that the gravitational back-reaction of an anisotropic string 
gas turns off at late times, due to the factor e*^ which appears in the equa- 

^To study the long-time behavior with an enhanced cross-section one should take ther- 
mal fluctuations into account, not only in setting the initial conditions as in H44|l . but also 
by adding an explicit noise term to the equations of motion. 
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tions of motion ((7j). One might be tempted to postulate a mechanism which 
stabihzes the dilaton, however this is problematic for reasons discussed in 
[23j : with pure Einstein gravity strings should freeze out, along the lines of 
our M-theory analysis 

Lest we appear too pessimistic, let us note some directions for future 
study which might invalidate our conclusions. 

• In this paper we have only studied decreasing dilaton solutions, whereas 
there is also a class of increasing dilaton solutions. If the value of the 
dilaton grew sufficiently large, the appropriate framework would be M- 
theory, and the results of our previous paper [2] would apply. However, 
it is possible that there is an intermediate time in which the dilaton is 
large enough for string annihilations to be effective, yet small enough 
for perturbative string theory to be relevant. 

• In this paper we made a number of simplifying assumptions. In par- 
ticular we assumed spatial homogeneity and only considered the radial 
moduli of the torus. A more complete analysis at the level of effective 
supergravity would be desirable; steps in this direction have been taken 
in [mini EH- 

The failure of the string gas scenario to naturally lead to three large spa- 
tial dimensions may be telling us one of three things. First, perhaps the string 
gas (or brane gas) framework is supplanted by other dynamics in the early 
universe, invalidating the approach we have been following. Second, per- 
haps the measure (jH^ does not reflect the true distribution of possible initial 
conditions of the universe. Third, perhaps three spatial dimensions is not 
favored. If one imagines that many "universes" are created, all with different 
initial conditions, then some sort of anthropic argument could be invoked. 
But many people, including ourselves, are uncomfortable with anthropic ar- 
guments until every other possibility has been explored. Consequently we 
intend to return to these cosmological issues as our understanding of string 
theory in the early universe improves. 
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Initial Number of Wrapped Dimensions 




Final Number of Wrapped Dimensions 




Figure 1: Average initial and final number of wrapped dimensions as a func- 
tion of the initial coupling and initial volume, starting with cp = —1 and 
evolved using the multiply-wound cross section. The volume is measured in 
units of {2TT\/a')^. 
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Initial Number of Wrapped Dimensions 




Final Number of Wrapped Dimensions 








Figure 2: Illustrates the dependence on the initial value of if. Same as Fig. 1 
except the simulations begin with (p = —0.15. 
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Figure 3: A histogram showing the distribution of the initial number of un- 
wrapped dimensions for universes that end up three dimensional. Extracted 
from the data set used to generate Fig. 1. 
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Figure 4: Histograms showing the distribution in the number of unwrapped 
dimensions at late times for various initial values of yj. Each histogram is 
based on 10^ simulations at an initial volume V = 4.0 x (2^T^/c^)^ and an 
initial ip = —1. 
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